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ABSTRACT 


Streamlined  shapes, both  two-dimensional  and  axi symmetric, are 
analytically  expressed  by  appropriate  polynomials  in  terms  of  in¬ 
dependent  parameters.  Permissible  ranges  of  the  independent  para¬ 
meters  are  examined  with  respect  to  selected  geometrical  constraints. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  sponsored  by  the  Naval  Ordnance 
Systems  Command  (Code  054131)  and  funded  under  UR  109-01-03. 

INTRODUCTION 

Streamlined  bodies  may  be  defined  as  those  bodies  with  negligible 
drag  due  to  separation  of  flow  on  the  tail.  Since  separation  is  a  vis¬ 
cous  flow  phenomenon,  its  occurrence  is  governed  by  the  Reynolds  number 
of  the  flow  past  the  body.  In  fact  at  extremely  low  Reynolds  numbers 
all  bodies  are  in  effect  streamlined  since  they  suffer  no  separation. 

However,  for  the  higher  range  of  Reynolds  numbers  of  practical  interest, 
it  has  been  experimentally  recognized  that  elongated  bodies,  that  is, 
bodies  with  rounded  noses  and  tapered  tails,  act  as  streamlined  shapes. 

In  this  respect  the  terms  elongated  bodies  and  streamlined  bodies  have 
often  been  used  synonymously  even  though  some  elongated  bodies  have 
proved  to  be  poor  streamlined  bodies. 

A  general  system  is  developed  for  analytically  defining  and  deter¬ 
mining  the  suitability  of  shapes  for  use  as  streamlined  bodies,  both  two- 


dimensional  and  axi symmetric .  Rounded,  pointed,  cusped, and  flat  faced 
ends,  with  and  without  parallel  middle  bodies,  in  any  combination  of 
forebody  and  afterbody  are  included.  Known  streamlined  shapes  may  be 
fitted  into  the  system  by  least-squares  fittings.  Polynomial  expres¬ 
sions  are  used  which  lend  t  emseives  most  readily  to  automatic  compu¬ 
tation. 

The  shapes  of  elongated  bodies  have  been  devised  by  potential  flow 

methods  such  as  source-sink  distributions  ^  in  order  to  obtain  knowledge 

of  the  flow  field  such  as  the  pressure  distribution  on  the  body.  However, 

2 

modern  methods  of  numerical  calculation  by  high-speed  computers  pro¬ 
vide  pressure  distributions  for  bodies  of  arbitrary  shapes  quite  readily 
so  this  aspect  is  no  longer  a  prime  consideration. 

The  devising  of  streamlined  shapes  has  had  Its  greatest  impetus 

from  the  development  of  aircraft  wings  for  two-dimensional  figures  and 

of  dirigible  bodies  for  axisymmetric  figures.  Analytic  expressions 

were  first  used  for  defining  shapes  from  which  pressure  distributions 

were  determined.  It  was  then  decided  to  start  with  whac  wei  considered 

suitable  pressure  distributors  and  determine  the  shapes  afterwards.  The 

3 

results  are  tables  of  offsets  as  given  by  the  NACA  airfoil  series  for 

1 4 

two-dimensional  figures  and  by  You'  Y  ana  Young  for  bodies  of  revolu¬ 
tion  with  rounded  noses  and  pointed  t;  Is. 

Since  both  the  pressure  distribution  and  the  axial  distribution  of  the 
offset  of  the  body  control  the  boundary-layer  development  leading  to 

^References  are  listed  on  page  83. 
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separation,  the  specification  of  pressure  distribution  alone  does  not 


ensure  the  best  streamlined  body  per  se.  The  approach  of  this  paper  is 
to  specify  the  geometrical  shape  of  families  ol  bodies  whose  hydrodynamic 
suitability  can  then  be  determined  by  further  analysis  such  as  that  of 
the  boundary-layer  development. 

A  notable  family  of  streamlined  bodies  of  revolution  was  developed 
by  Landweber  and  Gertler  ^  which  have  rounded  noses  and  rounded  tails. "Qua¬ 
dratic"  polynomials  are  used  in  whi»h  the  square  of  the  radial  offset 
is  a  polynomial  function  of  axial  distance.  What  is  significant  is  that 
independent  polynomials  are  related  to  the  independent  parameters  deter¬ 
mining  the  shape.  This  technique  was  presented  earlier  by  Admiral  David 
£ 

Taylor  in  his  studies  of  the  mathematical  description  of  ship  hulls. 

The  method  of  factorial  analysis  is  also  used  by  Landweber  and  Gertler 
for  determining  the  coefficients  of  the  independent  polynomials.  In 
addition  the  limiting  ranges  of  suitability  of  the  independent  para¬ 
meters  are  analyzed  by  an  envelope  technique  which  was  also  previously 
described  by  Admiral  Taylor.  The  "quadratic"  polynomial  developed 

from  delineations  of  dirigible  shapes  and  is  presented  in  a  primitive 
form  by  Lyon.  ^ 

The  "square  root"  polynomial  consists  of  an  ordinary  polynomial 

with  the  addition  of  a  square  root  term  to  accomodate  rounded  noses. 

It  was  used  In  the  specification  of  the  4  &  5-dfgit  series  of  airfoil 

89  1011 

sections  by  NACA  ’  and  by  Kwlk  ’  for  rudder  sections. 

The  method  of  this  report  follows  in  part  the  technique  of  Land¬ 
weber  and  Gertler  J  with  the  initial  difference  of  dividing  the  body 
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Into  a  fore  body  arid  afteL  body  a. id  then  normalizing  the  coordinate 
system.  The  result  is  that  there  are  only  two  independent  parameters 
remaining  to  consider  for  the  same  degree  of  generality.  Also  both, 
the  "quadratic"  polynomial  and  "square  root"  polynomial  are  subjected 
to  the  analysis  by  division  into  independent  polynomials;  the  factorial 
analysis,  where  allowable,  for  determining  the  coefficients  of  the  in¬ 
dependent  polynomials;  and  the  envelope  analysis  for  establishing  the 
range  of  values  of  the  independent  parameters  for  suitability  as  desir¬ 
able  shapes.  The  results  are  analytical  expressions  for  the  cases  of 
rounded,  pointed,  and  cusped  ends,  with  and  without  parallel  middle 
bodies,  in  any  combination  of  forebody  and  afterbody.  The  suitability 
requirements  are  presented  on  charts  giving  permissible  ranges  o£  the 
two  independent  parameters.  If  there  were  more  than  two  remaining 
independent  parameters,  such  a  simple  display  would  not  be  possible. 

The  special  case  of  the  flat  faced  nose  with  hydrodynamic  continuity 
is  finally  analyzed  by  a  "cubic"  polynomial. 
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GENERAL  ANALYSIS 


The  shapes  of  families  of  bodies  of  revolution  and  of  two-dimen¬ 
sional  symmetrl cal  bodies  may  be  stated  functionally  as 


Y  =  f  (X :  a£  ;  ) 


[I! 


where  Y 


X 


a 


1 


is  the  radius  of  the  body  of  revolution  or  the  offset  of  the 
two  dimensional  body, 

is  the  axial  distance  of  the  body  measured  from  the  nose, 
are  the  parameters  to  be  varied  which  specify  the  family  such 
as  length  L  ,  maximum  diameter  or  thickness  D  ,  etc.,  and 
are  the  boundary  conditions  or  restraints  which  give  desirable 
contours  such  as  closed  body  conditions  Y  =  0  at  X  =  0  and 
at  X  =  h  ,  etc. 


The  complexity  of  the  analytical  analysis  is  greatly  reduced  by  con¬ 
sidering  the  nose  and  tail  portions  separately  in  a  normalized  coordinate 
system  (>:,y):  y  =  0  at  x  =  0  and  y  =  1  at  x  =  1  .  For  the  curved 
body  the  split  into  forward  and  after  bodies  is  made  at  the  position  of 
maximum  radius  or  thickness.  For  the  parallel  middle  body  the  cuts  are 
made  at  both  ends  of  the  parallel  portion  of  the  body.  For  the  flat 
faced  body  an  additional  cut  is  made  at  the  edge  of  the  flat  face. 

For  the  completely  curved  body,  if  the  axial  distance  to  the  maximum 

section  is  X  ,  then  the  normalized  coordinates  become 
m 
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and 


2Y 


Y 

x^.  =  ^r~  for  the  forward  portion 

'  ni 


1  -  X 

x  ■=  - - —  for  the  after  portion 

a  L  -  X 


i?j 

[3] 

[M 


For  the  parallel  middle  body,  if  X  and  X  are  the  X-coordinat.es 

tl  t 

of  the  end  of  the  nose  and  the  beginning  of  the  tail  respectively,  then 
the  normalized  coordinates  become 


y 


2Y 

D 


[5] 


and 


x 


for  the  nose 


[bj 


L  -  X 

L  -  X. 


for  the  tail 


[7] 


For  the  flat  faced  body  the  normalized  coordinate  y  becomes 


y  = 


[8] 


where  D^.  is  the  diameter  or  width  of  flat  face. 

To  achieve  "hydrodynamic  continuity"  as  contrasted  with  mathematical 

continuity  it  is  only  necessary  that  the  position,  slope, and  curvature 

match  at  the  junction  of  the  forward  and  after  bodies.  Since  at  the 

dy 

junction  y  =  1  and  ^  =  0  ,  it  is  evident  that  the  position  and  slope 
requirements  are  always  met.  The  curvature  condition  however  remains  to 
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be  satisfied.  At  the  junction  the  curvature  is  given  by 
normalized  coordinates 


dX2 


X  .2 


'f  ‘(irhr) 


or 


h-Ci^) 


[10] 


where  k^.  =  curvature  of  forebody  at  x  =  1 


k^  =  curvature  of  afterbody  at  x  =  1 ,  and 
X 

m  =  —  ,  relative  axial  position  of  maximum  section. 


For  parallel  middle  bodies  it  is  obvious  that  k.  •=  k  =  0 

fa 


In  normalized  coordinates  the  contour  is  given  by 


y  =  f(x;at;fl.) 


i  =  1,  2, 

j  =  1,  2, 


[ii: 


where  a,  and  B.  are  now  defined  in  normalized  coordinates, 
i  J 

In  this  study  only  two  independent  parameters  and  are  to  be 

considered  for  simplicity  of  analysis. 

If  a  functional  form  like  that  of  a  polynomial  is  selected  as 


n=n 

y  =  Y  a  x11  =  P  (x) 
L>  n  n 

n=0 


[12] 


a  resolution  into  linearly  independent  polynomials  may  be  obtained  like 
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that  of  linearly  Independent  vectors  such  as 


>• -X-i  rn.i<,rt  "Iv-.j00  1131 

i  j 

or 

y  =  £  &1  Pn  1<x)  4  Q(x)  [14] 

l 

where 

Q(x)  =Y  0.  V  (x)  [13] 

1— i  J  n !  J 

j 

This  permits  the  effect  of  the  controllable  parameters  to  be 

obtained  Independently  of  each  other. 

The  Independent  polynomials  may  be  determined  by  substituting  con¬ 
ditions  and  /IT  into  the  general  polynomials  and  evaluating  the 

polynomial  coefficients  by  a  solution  of  the  resulting  simultaneous 
algebraic  equations  and  then  by  a  gathering  of  terms.  Another  method 
is  to  use  the  factorial  properties  of  polynomials  which  is  Illustrated 
in  the  specific  cases  to  follow. 

Not  all  variations  of  the  produce  desirable  shapes.  Conditions 

for  zero  values,  values  of  one,  maxima  or  minima,  and  inflection  points 
may  be  investigated.  For  example  the  condition  for  zero  values  of  y 

y(x;a^,a2>  =  0  ,  0  <  x  £  1  [16] 

may  be  studied  as  follows. 
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If  the  and  are  now  considered  as  variables  and  x  as  an 

adjustable  parameter,  a  line  may  be  defined  for  each  x  .  An  envelope 
to  these  lines  may  be  developed  which  represents  the  boundary  of  regions 
for  values  of  and  a0  with  and  without  an  additional  zero  value  of 

y  that  at  x  =  0  .  The  envelope  condition  is  given  by 


^  y(*;a1(V  -  0 


From  Equations  [16]  and  [17] 


“  ^(x) 


ci2  =  £  2  (x) 


A  plot  of  against  a2  for  the  range  of  values  of  x  , 

0  £  x  ^  1  ,  gives  the  envelope  curve. 
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"QUADRATIC”  POLYNOMIAL  REPRESENTATION 

GENERAL 

Tilt:  functional  relation 

n=n 

Y2  -  V  anxn  |20] 

n=0 

is  to  be  called  the  "quadratic"  polynomial  for  want  of  a  better  name.  It 
is  very  suitable  ior  describing  bodies  of  revolution  for  which  volume  is 
an  important  consideration  since  it  represents  the  axial  distribution  of 
the  cross-sectional  area.  It  has  the  additional  advantage  of  providing  a 
means  of  accommodating  the  analytical  description  of  bodies  with  rounded 
ends;  something  the  ordinary  polynomial  cannot  do. 

In  addition  to  bodies  with  rounded  ends  the  "quadratic"  polynomial 
may  be  applied  to  bodies  with  pointed  ends  and  cusped  ends.  Of  course 
since  the  representations  apply  equally  well  to  noses  and  tails,  any 
combination  can  be  formed  such  as  bodies  with  rounded  noses  and  pointed 
tails,  etc. 

Although  any  number  of  adjustable  parameters  may  be  used,  the 

analysis  is  to  be  limited  to  two  as  being  sufficient  for  describing 
suitable  figures.  Two  adjustable  parameters  for  each  partial  body  in 
the  normalized  coordinate  system  is  equivalent  to  six  adjustable  para¬ 
meters  for  the  whole  body  in  the  natural  coordinate  system. 
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ROUNDED  ENDS 


Since  there  are  five  conditions  in  all,  n  =  4  . 

The  and  ft are  substituted  into  the  polynomial ,  Differ¬ 

entiating  Equation  [20]  successively  with  respect  to  y  gives 
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[24] 


2y  *  (.flj  4 


2a7x 


3a3x 


4a(x3) 


dx 

dy 


and 


,  _  ,  2  .  3.  d^x 

(a  +  2a,x  4  3a  x  4  4a, x  )  — ^ 

“  14  dy~ 


4  (2a, 


[25] 


Since  a,  ^  0  ,  ~  “0  at  x  ■  0  .  This  automatically  provides  a 
1  o>  J  1 

rounded  end.  Then 


V 


2r 


The  other  substitutions  yield 


a 9)  2e  4  6a  +  12a.  =  -  2 k, 

2  2  3u  1 

fll>  ao  =  ° 

p2)  ao  +  ax  +  a3  4  a4  =  1 
(3^)  aj  4  2a2  4  3a3  4  4a^  =  0 

The  solution  of  Equations  [26]  by  determinants 

2 

a's  are  linear  functions  of  r  and  ,  Hence  y 

function  of  r  and  and  may  be  written  as 


[26] 


shows  that  the 
is  also  a  linear 


2 


y 


r  R(x)  4  k^K^Cx)  4  Q(x) 


[  27  J 


where  R(x)  ,  K^(x),  and  Q(x)  are  also  polynomials  of  the  fourth  degree 
in  x  .  It  is  possible  to  determine  R(x)  ,  K^(x),  and  Q(x)  by  first 
solving  for  the  a's  from  the  simultaneous  equations  and  then  regrouping 
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Another  method  is 


t emu  applicable  to  K(r.)  .  K,  (n),  and  ti<-v)  . 
developed  in  Reference  5  by  utilizing  the  factorial  properties  of 
polynomials  as  follows. 


It  is  evident  that  Che  relations  fot 


and  fi.  correspond  to 


V 

l  *2<»> 

7 

d"  2 

0..,) 

—  y  (i) 

dx 

7 

V 

X 

O  1 

O 

II 

O 

ft2) 

y“(D  =  1 

d  ■> 

183) 

Since  the  foregoing  apply  identically  to  r  and  kj  ,  it  is  further 
evident  that 


V 

R' 

(0) 

“  2  ,  KJ  (0)  =  Q'  (0)  -  0 

V 

V*» 

K1 

(D 

1! 

I 

to 

73 

l! 

O 

/“X 

>— 1 

II 

0 

V 

D 

x\ 

(0)  ■ 

=  Kx  (0)  =  Q  (0)  -  0 

Q 

(1)  = 

=  1  ,  R  (1)  =  Kj_  (D  =  0 

*3> 

R' 

(1) 

=  Kj  (1)  =  Q'  (1)  =  0 

2 

where  R'  *  ~  ,  R"  =  >  etc. 

dx  .2 

dx 

Evaluation  of  R(x) 


Since  R(0)  =  R(l)  =  R' (1)  =  R"(l)  =  0  and  R(x)  is  a  polynomial  of 
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the  fourth  degree,  R(x)  may  be  written  as 


Rlx)  =  ax  (x  -  I)"*  [30] 

Since  R1  (0)  «  2.  ,  a  =  -  2  . 

Then 


R('t)  =  -  2x  U  -  l)3 


[31] 


Evaluation  of  (x) 

A. 

Since  =  K? (0)  =  K^(l)  =  K^(l)  =  0  ,  K^(x)  may  be  written 

factorlaily  as 


Kt(x)  =  /?x2  (x  -  l)2  [32] 

Since  K'^  (1)  =  -  2  t  0  =  -  1  . 

Then 


K^(x)  =  -  x“  (x  -  l)2 


[33] 


Evaluation  of  Qfx) 

Since  Q'(0)  =  Q' (1)  =  Q"(l)  =  0  ,  Q'(x)  may  be  written  factorlaily 

as 

Q' (x)  =  yx  (x  -  l)2  [34] 
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Then  Integrating 


Q(x) 


2xJ 


3 


+ 


+  C 


With  Q ( 0 )  =  0  and  Q(l)  =  1  ,  C  =  0  and  y  =  12  .  Then 


Q(x)  =  x2(3x2  -  8x  +  6) 

For  rounded  ends  in  summary 

y2  =  r  R(x)  +  k^K^ (x)  +  Q(x) 

with 

R(x)  =  -  2x  (x  -  l)3 
(x)  =  -  x2  (x  -  l)2 
Q(x)  =  x2(3x2  -  8x  +  6) 

As  a  check  let  r  =  k^  =  1  which  are  the  conditions  for 

Then 


y 


2 


2 

x 


+  2x 


which  is  the  contour  of  a  sphere  with  center  at  x  =  1  . 

Graphs  of  R(x)  ,  K^(x),  and  Q(x)  are  given  in  Figure  1 


135] 

[36] 

[27] 

[31] 

[33] 

[36] 

sphere. 

[37] 
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Permissible  Ranges  of  Parameters  r  and  k( 

Not  all  combinations  of  r  and  kj  give  desirable  shapes.  Al¬ 
though  the  fourth  degree  polynomial  does  not  lend  itself  very  easily 
to  peculiar  shapes  for  positive  r  and  k^  ,  it  is  interesting  to  ana¬ 
lyze  the  possible  limitations  in  terms  of  simple  criteria: 

2 

1.  Zero  condition,  y  =0  for  0  s  >:  1  . 

2 

Negative  values  of  y  would  be  meaningless.^ 

2 

2.  Unity  condition.  y  =1  for  0  ^  x  £  1  . 

Bulges  above  y  =  1  are  undesirable. 

dy 

3.  Maximum  or  minimum  condition,  -r-  =  0  for  0  £  x  £  1  . 

dx 

No  other  maximum  or  minimum  is  to  be  permitted  than  at  x  -  1  . 

d2y 

4.  Inflection  point  condition.  — ~  ~  0  for  0  £  x  ^  1  . 

dx 

Inflection  points  are  undesirable  on  noses. 

Zero  Condition. 

y2  =  fCxjrjkp  =  0  ,  0  £  x  ^  1  [38] 

The  envelope  in  r  and  k^  with  x  as  the  variable  parameter  is 
given  by 

139) 

The  two  envelope  conditions.  Equations  [38]  and  [39],  provide  two 
simultaneous  equations  in  r  and  k^  which  are  solved  by  Cramer's  rule  to 
give  r(x)  and  kj(x): 
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r 


2 . 

*  l*  - 


[40] 


k 


1 


4x  +  6 


[41] 


The  envelope  curve  Is  shown  in  Figure  2.  Desirable  values  of  r  and 
k^  are  on  the  "inside  curved"  side  of  the  envelope  curve. 


Unity  Condition.  The  unity  condition  is  that 


y 


2 


f (x; r,k^)  =  1 


0  *  x  ^  1 


[42] 


The  envelope  in  r  and  k^  with  x  as  the  variable  parameter  is 
given  by 


(f  -  1)  =  0 


[43] 


The  two  envelope  condi tions , Equations  [42]  and  [43],  provide  two  simul¬ 
taneous  equations  in  r  and  which  are  solved  by  Cramer's  rule  to  give 


r 


1  + 


[44] 


k 


1 


[45] 
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Eliminating  x  gives  simply 


\  =  (2  -  r)2  [46] 

The  envelope  curve  is  shown  in  Figure  2.  Desirable  values  o£  r  and 
kj  ,  that  is,  without  bulges,  are  on  the  "inside  curved"  side  of  the 
envelope  curve. 

Maximum  or  Minimum  Condition.  The  maximum  or  minimum  condition  is 
given  by 


g  -  f  -  0  l«J 

The  envelope  curve  in  r  and  with  x  as  the  variable  parameter  Is 

given  by 

f"  =  0  [48] 

The  envelope  curve  is  shown  in  Figure  2.  A  better  understanding 
of  the  envelope  curve  is  developed  in  Figure  3.  Each  point  on  the 
envelope  curve  represents  a  tangent  giving  the  locus  of  values  of  r 
and  k^  which  provide  a  maximum  or  minimum  at  each  value  of  x  other 
than  the  maximum  at  x  =  1  which  prevails  at  all  times.  Two  such  loci 
are  represented.  Their  point  of  intersection  provides  a  value  of  r  and 
representing  maxima  or  minima  at  two  values  of  x  .  Evidently  from 
any  point  in  the  region  outside  the  envelope  curve,  two  tangents  may  be 
drawn  to  the  envelope  curve.  Thus  the  region  outside  the  envelope  curve 
represents  values  of  r  and  k^  giving  two  maxima  or  minima.  The 
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inside  the  envelope  curve  provides  no  maximum  or  minimum. 
Finally  there  is  only  one  maximum  or  minimum  specified  by  the  envelope 
curve  itself. 


The  two  envelope  conditions,  Lquations  [47]  and  [48],  provide  two  simul¬ 
taneous  equations  in  r  and  k^  which  arc  solved  by  Cramer's  rule  to  give  r(x)  and 


k^(x)  as 


r  = 


6x 


[49] 


6x  -  4x  +  1 


k  =  - 
K1  2 
1  6x  -  4x  4  1 


[50] 


Inflection  Point  Condition.  The  inflection  point  condition  is  given  by 


0 


For 


v2  =  fix) 


2f f"  -  f'2  =  0 


[51] 


and  the  envelope  condition 


f'"  -  0 


[52] 


The  two  conditions  provide  two  simultaneous  equations  in  r  and  in 

terms  of  x  .  Since  the  boundary  condition  leads  to  a  quadratic  relation. 


19 


Cramer's  rule  does  not  apply.  l'or  specified  values  of  x  the  two  simul¬ 
taneous  equations  may  be  solved  by  direct  substitution  of  one  equation 
Into  the  other.  The  results  arc  shown  in  Figures  2  and  4. 
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POINTED  ENDS 


The  adjustable  parameters 


are 


cij )  s  =  slope  at  x  =  0  ,  s  v  0 


\  dx  /y.=0 


[53] 


a^)  k^  =  curvature  at  x  =  1  ,  0 


k 


1 


dx2  ** 


[54] 


The  sign  of  the  equation  for  k^  is  chosen  so  as  to  require  to  be 

positive  for  desirable  shapes. 


The  boundary  conditions  /3^  are 


J3X)  x  =  0  ,  y  =  0 


fl9)  x  *  1  ,  y  =  1 
S3)  «  -  1  .  g  •  0 


[55] 


Since  the  "quadratic"  polynomial  automatically  gives  Infinite  slope  at 
x  =  0  ,  an  additional  condition  is  necessary  to  give  controlled  slopes 
at  x  =  0  .  Hence  the  degree  of  the  polynomial  becomes  five. 

For  a^) 
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(56] 


iiL 

dx 


a. 

1 


2a. x 
z 


+  . , 


na  x 

it 


n-1 


2y 


Since  y  “=  0  at  x-0,  ^  -  <*  unless  al  -  0  .  For  ^  -  0  , 

5^-  ig  indeterminate  at  x  “  0  .  Then  by  L  Hopital  s  Rule 
dx 


L 

x-,0 


±L 

dx 


(57] 


or 


z 


s 


=  a 


2 


then  requires  that 


a 


1 


0 


a 


2 


2 

s 


[58] 


The  other  substitutions  yield 


a2) 

2&2  +  + 

12a4  +  20a3  = 

V 

a  -  0 

0 

b2) 

a  +  a.  +  a 
o  1 

2  +  a3  +  a4  +  ‘ 

b3) 

al  +  ^a2  + 

3a„  +  4a,  +  5a 
3  4 

[59] 


y2  is  then  a  linear  function  of  s  and  k^  or 


y2  =  s2  S  (x)  +  (x)  +  Q(x) 


[60] 
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It  is  evident  that  the  relations  for 


and  Pj 


correspond  to 


s  d  2 

V  ^  y  'u)  " u 


~~~2  y  a) 

dx 

d2  7 

a2) 

0  i.  ,  m  v 

— n  y  0) 

dx“ 

¥ 

y2  (0)  =  o 

y2  0)  =  l 

*3> 

d  2  „* 
d^  y  (1) 

[61] 


Since  the  foregoing  apply  identically  to  s  and  ,  It  is  further 
evident  that 


aL)  S  *  (0)  =  K'(0)  =  Q’(0)  =  0 

S"(0)  =  2  ,  K” (0)  =  Q" (0)  =  0 
a  )  K"(l)  =-2  ,  S"(l)  -  Q" (1 )  =  0 

162] 

S(0)  =  K^O)  =  Q(0)  =  0 
fl2)  Q(l)  =  1  ,  S(l)  =  K^l)  =  0 
/33)  S'(l)  =  K’(l)  =  Q'(l)  -  0 


Evaluation  of  S(x) 

Since  S(0)  =  S ' (0)  =  S(l)  =  S'(l)  =  S"(l)  =  0  and  S(x)  is  a 
polynomial  of  the  fifth  degree,  S(x)  may  be  written  factorially  as 
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S  (x)  —  ax  (x  *  1) 

Since  S"(0)  »  2  ,  a  =  -  1  . 

Then 


S(x)  »  -  x2  (x  -  l)3 


lb  3] 


[64] 


Evaluation  of  (x) 

Since  Kl(0)  -  K[  (0)  =  K'^0)  -  K^l)  =  K[(l)  =0  ,  (x)  may  be 

written  factorially  as 

K^x)  -  <3x3  (x  -  l)2  [05] 

Since  K^'(l)  ■  -  2  ,  /3  =  -  1  - 

Then 

K^x)  =  -  x3  (x  -  l)2  [66] 


Evaluation  of  Q(x) 

Since  Q'(0)  =  Q"(0)  =  Q' (1)  =  Q"(l)  =  0  ,  Q(x)  may  be  written 
factorially  as 


Q'  (x) 


2  • 

Yx  (.x 


1)' 


e  r  -»  ' 

IV/  . 


Then  integrating 
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[68] 


With 

Then 


with 


Qw 


Q(0)  =  0  and  Q(l)  =1  ,  C  =  0  and  Y  =  30 


Q(x)  =  x3(fcx2  -  15x  +  10) 


[69] 


For  pointed  ends  in  summary 


y2  =  s“  S(x)  +  kj,Kj  (x)  +  Q(x)  [60] 


[64] 
[66] 
[69] 

Graphs  of  S(x)  ,  K^(x),  and  Q(x)  are  given  in  Figure  5  , 


S(x)  -  -  x2  (x  -  l)3 
Kl(x)  =  -  x3  (x  -  l)2 
Q(x)  =  x3  (6x2  -  15x  +  10) 
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Permissible  Ranges  of  Parameters 


and 


Four  conditions  are  to  be  applied:  the  zero  condition,  the  unity 
conditio. ,  the  maximum  condition, and  the  inflection  point  condition. 


Zero  Condition.  The  envelope  curve  is  specified  by 


and 


f (x;s2,k1)  =  0 

i  • 


0  ■£  x  ^  1 


The  results  are 


and 


x*  (3x  -  5) 
(x  -  1) 


3x'1"  -  lOx  +  10 
(x  -  l)2 


The  envelope  curve  is  plotted  in' Figure  6. 


Unity  Condition.  The  envelope  curve  is  specified  by 


[70] 


[71] 


[72] 


[73] 
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lOx  +  3 


(80] 


[81] 


=  30 (x  -  l)2 

10x2  -  lOx  +  3 


The  envelope  curve  Is  plotted  In  Figure  6. 


Inflection  Point  Condition.  The  enevelope  curve  is  specified  by 


2ff"  -  f'2  =  0 


[82] 


and 


£'"  =  0  .  [83] 

2 

The  variation  of  s  with  k^  is  obtained  numerically  by  direct  sub¬ 
stitution  in  solving  the  nonlinear  simultaneous  equations  of  the 
envelope  curve. 


For  x  =  C  an  indeterminate  condition  exists.  By  L'Hopital's  Rule 


[84] 


and 


3 (15s2  +  9k t  -  70) 
s 


The  boundary  curve  at  x  *  0  Is  then 


[85] 
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[86] 


or 


3s  -  k  +  10  =  0 


[87] 


The  point  of  tangency  at  x  =  0  Is  given  by  Equation  [87]  and 


(  ^  ‘ 0  ot  15s2  +  9ki  -  70  * 0  ■ 


[88] 


Then 


^  and  =  5  . 


[89] 


The  envelope  curve  is  shown  in  Figures  6  and  7. 
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CUSPliD  LNL'S 


The  adjustable  parameters 


arc 


cl.)  k  =  curvature  at  x  =  0  ,  k  '  0 
L  o  ’  o 


k 

o 


a9)  kj  =  curvature  at  x  =  1  ,  k^  "  0 


k 


1 


The  boundary  conditions 


H  are 
J 


«j_)  x  =  0  ,  y  =  0 

»2>  *  "  0  •  £  '  0 

x  =  1  ,  y  =  1 

R.  )  x  =  1  ,  ^  =  0 
4  dx 


[90] 


[91] 


[92] 


As  will  be  shown  the  cusped  end  requires  two  additional  conditions  which 
makes  n  =  7  . 


For  R  •  a  =  0 
1  o 

For  B2: 


*2  = 
dx 


a,  +  2a„x  +  3a„x  +  ...  +  na  x-1 
_1 _ 2 _ 3 _ n 


4^-  =  0  at  x  =  0  ,  a  =  0  . 
dx  1 


[93] 
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.  2/1,  4  j.l  .  .  .  4  |l(|  X 

Jy  2  i  i> 

Then  ,  a  — - — — — - 

<iX  7  k 

2  <"2  *  V  *  ...  4  )4 


(  £  I-„ 


F"'  (  S  i-o  ■  "  •  •,  ■ » 


For  a.  : 


,2  6a, x  4-  12a, 4  ...  ■*•  n(n  -  l)a  x°  ^ 

<1  y  _ 3 _ 4  n _ 


,  , .  2  .  n-2  2  .  2  .  3  n-1,2 

6a_x  4  12a. x  ♦  ...  ♦  nCn  -  l)a  x  (3a, x  *  4a,x  ■*  4a. x  ♦  . . .  ♦  na  x  J 

34  v/n  v  3  4  4  n 


.  .  3  4  n .  k 

2(a.x  +  a  x  +  ...  +a  x  )  * 

3  4  n 


.  .  3  A  n.3/2 

A^X  +  fl^X  -»...+  B^x  ) 


2[6a,x  +  12a, x2  +  ...  ♦  n(n  -  l)a  xn"2J(a,x3  ♦  a.x4  ♦  ...  ♦  a  xn)  -  (3a, x2  ■*  4a, x3  ♦  ...  +  na  x11  ) 

3  4  n  3  4  n  3  4  n 


,  ,  3  4  ^  n.3/2 

4  («3x  +  a^x  +  .  .  .  -f  aRX  ) 


2[6a  ♦  12a.x  +  ...  +n(n  -  l)a  x°  ^(a.  4  ax  4  ...  +  a  x°  -  (3a  4  4a  x  4  ,  .  .  +  na  x°  ^ ) 

^34  n  3  4  n  34  n 


.  1/3  .  A/3  .  n-8/5,3/2 

V«3*  +  *4X  T  •••  +  anx  > 


At  x  =  0  , 


Let  «3  ■  0 


2[l2a,  +  20a  x  +  .  ,  .  +  n(n  -  lja  xn  **](a  +  a  x  +  —  +  a  xn  *)  -  (Aa,  +  Sa,x  +  ...  tm  i"  3)2 

*-4  5  n  4  j  n  4  o  n 


4(a.  .  „  n-4.3/2 

4  +  a5x  4-  . . .+  aRX  ) 


ii 


At  x  «  0 


£4  -  2a,  ' 
dxZ  4 


1102) 


4a,  =  k 
4  o 


[103] 


The  other  substitutions  yield 


a„-.  12a.  +  20ac  +  30at  =  -  2k, 

2  4  5  6  1 


«  •  a  +  a  +  a  =  1 
3  4  5  6 


[104] 


jS,  :  4a  +  5a  +  6a  =  0 
4  4  5  6 


2  2 
y  is  then  a  linear  function  of  k  and  k, 

o  1 


y2  =  kQ2  Kq(x)  +  kj  Kj (x)  +  Q(x) 


[105] 


It  is  evident  that  the  relations  for  and  correspond  to 


V  4  *2«»  ■ 0 

dx 


d*  2  2 

—  y  =  6k0 

dx 


j2 

aj  y  (D  «  -  2k 

*  dx  1 


PJ  y  (0)  «  o 


»2>  £  *2«»  ‘  0 


[106] 


d2  2 

—  yz  (0)  =  0 

dx^ 
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it  is 


i V  y"U)  -  1 

<V  ^  (1)  ■  0 

Since  the  foregoing  apply  identically  to  kQ  and  , 

further  evident  that 

1(1  III  III 

a.)  K  (0)  =  K.  (0)  =  Q  (0)  -  0 

1  o  1 

KoiV(0)  =  6  ;  KxiV  (0)  =  QiV  (0)  =  0 

a.,)  KV  (1)  =  -  2  ;  K"  (1)  =  Q"  U)  *  0 

2  1  o 

S1)  Ko(0)  =  K^O)  =  Q(0)  =  0 

[107] 

P2)  Kq'  (0)  =  K1'  (0)  =  Q1  (0)  =  0 

K"  (0)  =  KV  (0)  =  Q"  (0)  =  0 
o  i 

fi3)  Q ( 1 )  =  1  ;  Ko(l)  =  K^l)  =  0 

0.)  K  '  (1)  =  K.  '  (1)  =  Q  (1)  =  0 
4  o  i 

Evaluation  of  K  (x) 
o 

Since  K  (0)  =  K  '  (0)  =  K  "  (0)  =  K  (0)  =  K  (1)  =  K  '  (1)  = 
o  o  o  °  o  o 

K  "  (1)  =  0  and  K  (x)  is  a  polynomial  of  the  seventh  degree,  K  (x) 
o  o  o 

may  be  written  factorially  as 

Ko(x)  =  ax4  (x  -  l)3  [108] 

Since  K  1V  (0)  =  6  ,  a  =  -  1/4. 
o 
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Then 


K0(x)  “  -  (*  -  l)3  [109] 

Evaluation  of  K^(x) 

Since  Kl(0)  -  ^'(O)  =  ^"(0)  »  ^*"(0)  -  K  iv  (0)  -  K^l)  = 

K|'(l)  ■  0  ,  (x)  may  be  written  factorially  as 

K^x)  ■=  flx5  (x  -  l)2  [110] 

Since  K^'fl)  -  -  2  ,  0  =  -  1  . 

Then 

K^x)  =  -  x5  (x  -  l)2  [111] 


Evaluation  of  Q(x) 

Since  Q'(0)  =  Q" (0)  =  Q'"  (0)  =  Q1V  (0)  =  Q' (1)  =  Q"(l)  =  0  , 

Q1 (x)  may  be  written  factorially  as 

Q' (x)  «  yx4  (x  -  l)2  [112] 


Then  integrating 


7  6  5 

Q(x)  =y(^r-^_  +  2L.)+c 


[113] 


With  Q(0)  =  0  and  Q(l)  =  1  ,  C  =  0  and  y  =  105  . 
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f 


Q(x)  =  x5(15x2  -  33x  +  21) 


For  cusped  ends  In  summary 


y“  *  ko2  Kq(x)  +  \  K^x)  +  Q(x) 


K0(x)  =  -  ix4  (x  -  l)3 


K1(x)  =  -  x5  (x  -  l)2 


Q(x)  =  x5(15x2  -  35x  +  21) 


Graphs  of  Kq(x)  ,  K^x),  and  Q(x)  are  given  in  Figure  8. 


Permissible  Ranges  of  Parameters  k  and  k, 

o  1 


As  before  four  conditions  ore  to  be  applied:  the  zero  condition, 
the  unity  condition,  the  maximum  and  minimum  condition, and  the  inflec¬ 
tion  point  condition. 


Zero  Condition.  The  envelope  curve  is  specified  by 


f  =  0 


and 


f'  =  0 


[116] 


The  results  are 


W  2  =  4x2(5x  -  7) 
°  (x  -  l)3 


[117] 


and 


kl  = 


lGx 


28x  +  21 


(x  -  1)' 


[118] 


The  envelope  curve  Is  shown  In  Figure  9. 


Unity  Condition.  The  envelope  curve  is  specified  by 


f  -  1  =  0  [119] 

and 

f'  =  0  [120] 


The  results  are 


4(5x^  +  8x3  +  9x2  +  8x  +  5) 

4 

x 


[121] 


and 


k 


1 


10x5  -  8xA  -  5x3  -  2x2  +  4 
5 

X 


[122] 
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Hie  envelope  curve  is  shown  in  Figure  9, 


Maximum  or  Minimum  Condition.  The  envelope  curve  is  specified  by 


f  -  0  [123] 

and 

f"  =  0  [124] 


The  results  are 


420x2 

2ix2  -  28x  +  10 


[125] 


and 


210 (x  -  l)2 _ 

21x2  -  28x  +  10 


[126] 


The  envelope  curve  is  plotted  in  Figure  9. 


Inflection  Point  Condition,  The  envelope  curve  is  specified  by 


2ff" 


1127] 


57 


and 


£"'  =  0  [128] 

2 

The  variation  of  k  with  k,  is  ootained  numerically  by  direct 
o  1  J  1 

substitution  In  solving  the  nonlinear  simultaneous  equations  oi  the 
envelope  curve. 

The  results  are  shown  in  F.  ^ora  10.  The  various  regions  giving  the 
number  of  inflection  points  are  delineated  in  the  same  figure. 
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The  signs  ot  the  equations  for  r  and  k^  are  chosen  so  as  to  require 
r  and  to  be  positive  for  desirable  shapes. 


The  boundary  conditions 


are 


8y)  x  =  0 
B2)  x  =  1 
8^)  x  »  1 


y  =  0 


[132] 


Since  there  are  five  conditions  in  all,  n  =  3  . 

The  and  8 are  substituted  into  the  polynomial.  For  : 

Differentiating  Equation  [129]  with  respect  to  y  gives 


1 


(^a^x 


-h 


ai  + 


2a  2x  + 


3a3x2) 


dx 

dy 


[133] 


or 


dx 

dy 


Js 


\a,  +  a  x^  +  2a  x3/2  +  3a  x5/2 
1  2  3 


[134] 


dx 

At  x  =  0  .  —  =  0  which  ensures  a  rounded  end, 
dy 

[133]  with  respect  to  y  gives 


Differentiating  Equation 


/  t  "3/2 
(-^a^x  + 


2a  2  + 


6a3x) ( ~ )2  +  (%a^x*^ 


a^  +  2a2x  + 


3a3x2) 


dy 


[135] 


or 


d2x 

-7 


+  a^x^  +  2a2x3/2  +  3a3x5/2)  +  a^  +  2a2x3/'2  +  3a3x5/2)2 


^  -  2a2x3/2  -  6a3x5/2 


[136] 
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At  x  =  0 


,  l 
r  =  ^ 

or 

a^  =  'J  2r 

a^)  -ia^  +  2a?  +  6a3  =  -  k^ 

*1>  ao  =  0 

[137] 

S2}  %  +  al  +  a2  +  a3  =  1 

P;)  +  a1  +  2a2  +  3a3  =  0 

The.  presence  of  the  square  root  term  prevents  the  use  of  the  factorial 
analysis. 

The  solution  as  simultaneous  equations  in  the  a's  produces 


y  -•  v'  2r  R(x)  +  k^  (x)  +  Q(x) 

with 

R(x)  =  V~x  -  |  (3x2  -  lOx  +  15) 
K^(x)  =  -  4x  (x  -  l)2 
Q(x)  =  x  (x2  -  3x  +  3) 


[138] 

[139] 

[140] 

[141] 


Graphs  of  R(x)  ,  (x) f  and  Q(x)  are  shown  in  Figure  11. 
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Permissible  Ranges  of  Parameters  r  and 

For  the  same  reasons  as  the  "quadratic"  polynomial  the  conditions 
to  be  considered  are  the  zero  condition,  the  unity  condition,  the 
maximum  or  minimum  condition, and  the  inflection  point  condition. 

Zero  Condition.  The  envelope  conditions  are 


y  =  g(x;r.k^)  =0  0  ^  x  ^  1 


[142] 


and 


g-  =  0 


[143] 


Ihe  results  are 


V  2r  = 


2 (x  -  3) 


,  -1/2  -3/2 

5x  -  x  +  x 


[144] 


and 


kl  = 


-1/2  1/2  -3/2  2 

36x  +  20x  +  12x  +  x  -  12x  +  15 


2  (x  -  1) (5x 


-1/2 


-3/2  . 

X  +  X 


5) 


Ihe  results  are  plotted  in  Figure  12. 


[145] 


Unity  Condition.  The  envelope  curve  is  specified  by 
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42 


(153] 


f 


6 (x  -  1)  (  -  1  -t-  5x  4  x5/2  -  5x3/2) 

1  -  1 2x  +  15x2  +  12x7/2  -  36x5/2  +  20x3/2 


The  envelope  curve  Is  plotted  In  Figure  12, 


Inflection  Foint  Condition.  The  envelope  conditions  are 


g"  =  0  (154] 

and 

g'"  =  0  [155] 


The  results  are 


J~2x 


16x 


5/; 


Sx3/ 2  -  5x  +  2 


[156] 


and 


k 


1 


2(2x5/2  -  5x  +  3) 
5/2 

8x  -  5x  +  2 


[157] 


The  envelope  curves  are  plotted  in  Figures  12  and  13. 


POINTED  ENDS 


The  ordinary  polynomial  is  utilized 


y 


r 

a  x 
n 


[158] 
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The  adjustable  parameters  a,  arc 


a^) s  =  slope  at  x  =  0  ,  s  ^  0 


\  dx  /x= 0 


[159] 


o^)  k^  =  curvature  at  x  =  1  ,  &  0 


[160] 


The  sign  of  the  equation  for  k.  is  chosen  so  as  to  require  to  be 

positive  for  desirable  shapes. 

The  boundary  conditions  fi .  are 


«1> 

x  =  0 

»  y  =  0 

V 

X  =  1 

.  y  =  i 

[161] 

fl3) 

X  =  1 

1! 

O 

Since  there  are  five  conditions  in  all,  n  =  4  . 

Substitution  of  and  into  the  polynomial  produces 

al>  a^  =  s 

a„)  2a  +  6a  +  12a  =  -  k 

2  2  j  4  1 

*!>  ao  =  0  [162] 
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n2)  ao  +  +  a2  +  a3  +  a4  -  l 


By)  ^  +  2a £  +  3a3  +  4a^  -  0 


y  Is  a  linear  function  of  s  and  or 

y(x)  =  s  S(x)  +  kj  K^Ox)  +  Q(x) 


The  relations  for  and  correspond  to 


V 

y '  (0) 

=  s 

a2) 

y"  (1) 

=  “  ki 

V 

y  (0) 

=  0 

V 

yd) 

=  l 

V 

y'd) 

-  0 

evident 

that 

V 

<•  1  /rt\ 

5  W  = 

1  ;  KJ 

(0)  =  Q'  (0)  = 

a2) 

II 

i- 4 

a 

-  1  ; 

S" (1)  =  Q"(l) 

*L> 

S(0)  = 

K(0)  = 

Q(0)  =  0 

02> 

Q(  1)  = 

1  ;  S(l)  =  K^l)  =  0 

V 

S' (1)  = 

K((l) 

=  Q'(l)  =  0 

Since  S(0)  =  S(l)  =  S'(l)  =  S"(l)  =  0  ,  S(x)  =  ax(x  -  1) 
Since  S'  (0)  =  1  ,  a  =  -  1  and 

S (x)  =  -  x(x  -  l)3 


[163] 


[164] 


[165] 


[166] 
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Since  K^O)  -  Kj  (0)  -  K^'l)  =  K[(l)  =  0  ,  K^x)  =  /9x2(x  -  l)2  . 
Since  K'j’(l)  =  -  1  ,  jS  =  \  and 

K^x)  =  -  *.x2<x  -  l)2  [167] 

Since  Q'(0)  =  Q' (1)  =  Q"(l)  =  0  ,  Q' (x)  =  Yx(x  -  l)2  . 

Since  Q(0)  =  0  and  Q(l)  =1  ,  y  =  7  and  then 

Q(x)  =  x2(3x2  -  8x  +  6)  [168] 

In  summary  for  pointed  ends 


y  =  s  S(x)  +  kL  K1(x)  +  Q(x) 

[163] 

S(x)  =  -  x(x  -  I)2 

[166] 

Kl(x)  =  -  ^x2(x  -  l)2 

[167] 

Q(x)  =  x2(3x2  -  8x  +  6) 

[168] 

The  polynomials  are  plotted  in  Figure  14. 

Permissible  Ranges  of  Parameters  s  and 

Four  conditions  are  to  be  considered:  the  zero  condition,  the 
unity  condition,  the  maximum  or  minimum  condition,  and  the  inflection 
point  condition. 
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Zero  Condition.  The  envelope  curve  is  specified  by 


g  =  0 

an  i 

g'  =  0 


The  results  are 


The  envelope  curve  is  plotted  in  Figure  15. 


Unity  Condition.  The  envelope  curve  is  specified  by 


(169] 


[170] 


[171] 


[172] 


g  -  1  =  0  [173] 

and 

g’  =  0  [174] 
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Tli e  results  arc 


s  =  2(1  +  h  [175] 

and 

\  =  2(1  -  V  [176] 


Ihe  envelope  curve  is  plotted  in  Figure  15. 


Maximum  or  Minimum  Condition.  The  envelope  curve  is  specified  by 


S’  =  0  '177] 

and 

g"  =  0  [178] 


The  results  are 


s 


12x 


6x 


4x  +  1 


[174] 


and 


k 


1 


12(x  -  l)2 

2 

6x”  -  4x  +  1 


The  envelope  curve  is  plotted  in  Figure  15. 


[180] 
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Inflection  Point  Condition,  The  envelope  curve  is  specified  by 


[181] 


and 


8'"  =  0 


[182] 


The  results  are 


4(3x2  -  3x  4  1) 
6x“  -  8x  +  3 


1183] 


and 


6x2  -  8x  +  3 


[184] 


For  x  =  0  the  boundary  line  given  by  g"  =  0  is 


6s  +  kl  -  12  =  0  . 


[135] 


The  envelope  curves  are  shown  in  Figures  15  and  16. 


CUSPE1)  ENDS 


The  ordinary  polynomial  is  utilised 


SO 


n=n 

V  . 

Li  n 
n-0 


Tlic  adjustable  parameters  c.^  arc 


a,)  k  =  curvature  at  x  =  0  ,  k  '•  0 
i  o  o 


f  d  y  \ 
"°  “  -  dx2 


o.^)  k^  =  curvature  at  x  =  1  ,  kj  0 


The  boundary  conditions  0  are 


01)  x  =  0  ,  y  =  0 


Aj)  x  =  1  ,  y  =  1 

V  - ' 1  •  £  - 0 


Since  tht^  ■re  r.ix  conditions  :r.  all,  n  ---  3  . 


Substitution  of  and  0  into  the  polynomial  produces 


V  2a2  =  ko 

2a0  4  6a^  +  12a  4  20a  =  -  k 

2  2  3  4  5  1 


V  ao  -  0 


51 


Since  K  (0)  =  K  '(0)  =  K  (1)  “  K  '(1)  «  K  \'U  ■=  0  , 

0  0  o  o  o 

K  (x)  =  cxx"(x  -  l)-'  1 1 94 J 

o 

Since  K  " (0)  =1  ,  a  =  \  and 
o 

K  (x)  =  ijx""  (x  -  l)3  [155] 

o 

Since  K,  (0)  =  Kj'(O)  =  K  "■<())  =  1^(1)  =  Kj  '  (1)  =0  , 

K^(x)  =  )9x3  (x  -  l)2  [156] 

Since  K^'d)  «  -  1  ,  P  -  -  ^  and 

Kl(x)  =  -  ^xJ (x  -  l)2  [157] 

Since  Q(0)  -  Q"(0)  =  Q'd)  =  Q"(l)  =  0  , 

k 

Q' (x)  =  \x2(x  -  l)2  [198] 

Since  <j(U)  =  0  acd  Q(lt  =-1  ,  V  =  30  and 

Ci(x)  =  x3(6x2  -  15x  +  10)  .  [199] 

In  summary 
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y  -  k 

1  o 

K  (x)  4  k.  K. (x)  ^  Q(x) 
o  11 

[191] 

k  (x) 

0 

.  2.  . .  3 

*  XX  lx  -  i) 

[19b] 

K^x) 

3  2 

=  -  1.x  (x  -  1) 

[197] 

Q(x)  = 

x3(6x2  -  15x  +  10) 

[199] 

These  are  plotted  in  Figure  17. 


Permissible  Ranges  of  Parameters  kQ  and 

Four  conditions  are  to  be  investigated:  the  zero  condition,  the 
unity  condition,  the  maximum  or  minimum  condition,  and  the  inflection 
point  condition. 

Zero  Condition.  The  envelope  curve  is  specified  by 


8  =  0 


[200] 


and 


8'  =  0 


[201] 


The  results  are 


k 

o 


2 x  (3x  -  5) 
(x  -  l)3 


[202] 


and 
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|  203] 


k 


1 


v  x 


10x  -i  iv; 
2 

‘ )  “ 


Tlu;  envelope  enrve  Is  plotted  In  1 1 pure  18. 


Unity  Condition.  The  envelope  curve  is  specified  by 


8  -  1  =  0 


and 


8 


0 


The  results  are 


k 

o 


2(3x2  +  4y  +  3) 

2 

x“ 


and 


k 


1 


2  (x  -  1)  "  (3x  4  2) 
2 


x 


The  envelope  curve  is  shown  ir.  Figure  18. 


Maximum  or  Minimum  Condition.  The  envelope  curve  is 


[204] 


|205] 


[206] 


[207  ] 


specified  by 
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s’  -  C  [208] 

and 

g"  =0  [ 2C9  ] 


The  results  are 


k  = 
o 


(.'Ox 


10x“  -  lOx  4  3 


[210] 


and 


60  (x  -  l)*- 
10x2  -  lOx  4  3 


The  envelope  curve  is  plotted  in  Figure  18. 


[211] 


Inflection  Point  Condition.  The  envelope  curve  is  specified  by 


8 


M 


[?12] 


and 


8 


111 


0 


[213] 


The  results  are 


k  = 
o 


20x2(6x2  -  8x  4-  3) 
~Z  "  .  „  3  2 


-  20x  4-  4Cx  -28x  4-  8x  -  1 


[2141 


56 


“et  BK"  I'lHYV'MIA!  Ktd’KIlShN  ’  A  i  l«'N 
I'OR  !'1,A  i  !AOH>  NOSkS 


To  adiicvt  "hy Ji odynaml c  continuity"  nt  the  edge  ol  the  ll.it 
laced  nose,  ?.ctj  curvature  Is  req.ilrcJ  In  addition  to  Infinite  slope. 

A  "cuMe"  polynomial  achieves  this,  namely, 

3  .’3 

y  c  P(x)  =  a  -I  a,x  H  a.,x‘'  a  a  .x  -i  .  .  .  l?lo] 

o  1  d  o 


Since 


dx 

dy 


0 

_ 3v^ _ 

aL  *  2a2X  '  3a3"'2  4 


11117] 


and  y  »  0  sit  x  *=  0 


C  <Iy  ^X~0  ~  °  (infinite  slope) 


1218] 


Also  since 


6y 


(2a2  a  6a^x  -t  . 

o 

+  Za^x  -t  3  a  ^  x 


.  •  1 


r&i 

\  dv  J 


-t  .  .  . 


[219] 


and  7  0 


(zno  curvature) 


[220] 


There  Is  one  adjustable  parameter 
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k,  «  ruivatuu  at 


0 


1 

3l 


K  »  I  ,  k. 


[221] 


Hit  bound.-,  ry  coud  i  l  Ion*.  /Sj  arc 


x  “  0  ,  y  «  0 
/U  x  ■=  1  ,  y  «  1 

i. 

V  *  “  1  >  £  "  ° 


i22J 


flute  i litre  are  lour  conditions  in  oil,  the  polynomial  In  x  is 


cubic. 


!  'i 

arc 

subst i tutcJ  Ini  o 

al 1 

2a2 

*  °;i3 

"  '  3kl 

'V 

o  * 
o 

0 

/U 

a  ~t 
o 

*1  4 

a  „  h  a  „  *-  1 
2  3 

2a.,  -t 

3a3  =  ° 

I  -23  J 


The  form  of  the  "cubic"  polynomial  Is  then 


>•  ■=  kj  kj(x)  -t  Qtx) 


12  24] 


The  relations  for  a.  and  fi .  eoriespond  to 
*  i 


V 


C)  =  *  3k, 


dx* 
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[  225] 


V 

y  (0)  =  0 

y  3  ( 1 )  =  1 

V 

evi Jcnt 

1  hat 

Kyd)  =  -  3  ,  Q"(l)  =  0 

ft )  K.  (0)  =  Q(0)  =  0 

1  1  [22b] 

ft,)  Q  ( 1 )  »  1  ,  Kx  (1)  =  0 

ft3)  (1)  -  Q’  0)  =  0 

Evaluation  of  K^(x) 

Since  1^(0)  =  Kjd)  =  K^l)  -  0 

Kl(x)  =  ax(x  -  l)2  [227] 

Since  K'^d)  -  -  3  ,  a  =  -  |  •  Then 

Kl(x)  -  -  |  x(x  -  l)2  [228] 

Evaluation  of  Q(x) 

Since  Q"(l)  =  Q ’ (3 )  =  0 


Q'  =  0(x  -  l)2 


[229] 
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=  c 


Since  g(0)  =  0  and 


Q(1 )  =1  ,  £  =  3  and  C  =  1 


Q  =  (x  -  1)  i  1 


independent  polynomials  are  plotted  in  Figure  20. 


Permissible  Ranges  of  kj 


Zero  Condition. 


,3  *  h(x,k)  "  0  ,  0  x  •=  1 


The  boundary  condition  h  ®  0  gives 


k,  - 

1  3x(x  •  1)' 


For  x  =  0  ,  k 


is  indeterminate  and  by  l.'Hopital  s  Rule 


k  =  ° 
kl 


From  the  plot  of  ^  against  x  in  Figure  21  it  is  evident  that  the 


permissible  range  is 
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2 


[235] 


k 


1 


Unity  Condition. 


h(x,k)  s.  1  ,  0  •.  x  1 


[236] 


The  boundary  condition  is  then 


h  =  1 


1237  ] 


Then 


k 


1 


[238] 


From  the  plot,  of  kj  against  x  in  Figure  21  it  is  evident  that  the 
permissible  range  is 


k 


1 


;>  0 


[239] 


Maximum  or  Minimum  Condition. 


4^  >  0  ,  0  ^  x  <  1 

dx 


[240] 


except  at  x  =  1  .  Then 


h' 


x  -  1 


<  0 


[241] 
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Tlic  boundary  condition  it  then 


1. ' 


x  -  1 


0 


[242] 


which  produces 


k 


1 


2(x  -  1) 
3x  -  1 


[243] 


From  the  plot  oi  against  x  in  'figure  21  it  is  evident  that  the 

permissible  range  is 


0 


|244] 


Inflection  Point  Condition. 


0  <  x  •=  1 


[245] 


or 


3hh"  -  2h'2  <  0 


9h 


5/3 


[246] 


For  h  0  (positive  condition) 


3hh"  -  2h’2  =  C 


[247] 


which  reduces  to 
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(x  -  l)2k12  +  2 (3x  -  2)k  =  4 (x  -  1) 


3x  +  x-/  4x  -  3 
(x  ■■  l)2 


The  perrr.i edible  range  is  0  and  with  the  zero  condition  of  ^  2 

the  combined  permissible  range  is 


0  i  ^  2 


as  shown  in  Figure  21. 
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LEAST  SQUARES  FIT 


general 

Known  shapes  may  be  approximately  expressed  In  "quadratic" 
and  "square  root"  polynomials  by  a  least  squares  fit.  If 


z  =  y  for  "quadratic  1  polynomials 


or 


z  y  for  "square  root"  polynomials 


then  in  general 


z  =  d  D(x)  +  (x)  +  Q(x) 


[251] 


where  d  and  D(x)  refer  to  the  appropriate  type  of  figure.  For 
example  for  a  rounded  nose  shape  d  =  r  and  D(x)  =  R(x)  .  Since 
the  whole  body  is  to  be  fitted,  both  fore  and  after  bodies  are  to  be 
considered  which  would  mean  determining  four  coefficients  d^  ,  , 

d  ,  and  k  by  a  least  squares  fit.  The  subscript  f  refers  to 

the  forebody  and  a  to  the  afterbody.  However  since  k^  and 

k,  are  related  bv 

la 


[252] 
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BODIES  OF  REVOLUTION 

In  the  normalized  coordinates  of  this  paper  all  ellipsoids 
are  transformed  into  spheres  so  that  r  =  k^  =  1  which  is  plotted 
In  Figure  4  for  bodies  with  rounded  ends. 

An  elilpsoid-like  body  was  developed  by  Munzner  and  Relchardt^*- 
which  has  an  almost  constant  pressure  distribution.  The  Relchardt  body 
is  expressed  as 
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1  O  A 

(*-rO  <(0  " 


[257] 


whore  a  and  b  arc  the  semlaxcs  and  ir  normalized  coordinates  as 


■>  o  A 

(x  -  1)  +  y  ’  =  1 


[258] 


For  a  least  squares  fit 


Zj  =  [x(2  -  x) ] 


1 

1.2 


[255] 


A  least  squares  fit  by  rounded  end  "quadratic"  polynomials  gives 


r  =  1.344 


and 


k  =  1.085 


which  is  plotted  in  Figure  4. 

Other  known  bodies  of  revolution  may  be  fitted  and  plotted  in 
Figure  4  for  comparison. 

TWO-DIMENSIONAL  SHAPES 

In  normalized  coordinates  all  ellipses  have  been  transformed 
into  a  circle 
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2.08.’ 


and 


-  2.437 

which  arc  plotted  In  V inures  13  and  10. 

1  3 

The  Navy  Standard  Strut  was  superseded  by  the  Id'll  shape 
which  has  a  rounded  tall.  A  least  squares  fit  yields  for  the  nose 
port i on 


•J~~ r  =  1.41A 


and 


kf  =  0.752 


and  for  the  tall  portion 


J~lr  =  0.389 


and 


kt  =  1.124 


which  are  plotted  in  Finn re  13. 
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1’RISMAI'U:  C0I.1TU  1  UN  T 


The  fullness  ol  a  shape  Is  given  hy  the  prismatic  coefllctenl 
w,.,cll  thc  ratio  oi  the-  volume  ol  the  body  to  the  volume  of  a  prism 
having  the  maximum  o  loss  -  see  t  Lon  area  and  the  lcupth  o)  ihe  lods 
lor  a  body  oi  revolution  the  prism  Is  represented  by  a  cylinder  ol 
ma x Lmum  diameter  aiul  ol  tdie  Vn>d^ 

The  prismatic  cod' Tic  lent  of  the  semibodies  ol  this  report  is 
given  by 

c  -  f'  >h«  1^1 

P  Jo 

tlic  prismatic  coefficient  of  the  whole  body  is  given  by 

C  =•  r-  C  ■»  0  -  m)  C  I2o3J 

1>  l'f  ta 

where  the  subscripts  f  and  a  refer  to  the  forebody  and  afterbody 
respectively. 

The  utility  of  the  quadratic  polynomial  in  computing  the  pris¬ 
matic  coefficient  is  evident  from 


„i  rl  r1 

Cp  =  dI°  D^dx  +  kl]o  Ki^dx  +  jo  ^dx 


For  the  rounded  end 


[  265] 
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iiiiu  Lhc  pointed  uiu 


0 


r 


o 

s‘ 

00 


[206] 


and  the  cnsped  end 


C 

V 


1120 


[  207] 


Lines  of  constant  <J 

r 


arc  ploLLed  In  l'ijjui'os  4,  7,  and  10. 


i 


Figure  3  -  "Quadratic"  Polynomial:  Rounded  Hint  -  Delineation  of 
Regions  by  Lnvelopc  Curve  for  Maximum  and 
Minimum  Condition 


Figure  4  -  "Quadratic"  Polynomial:  Rounded  Fnd  -  Inflection  Point  Condition 


Figure  7  -  "Quadratic"  Polynomial:  Pointed  lind  -  Inflection 

Point  Condition 


Figure  8  -  "Quadratic"  Polynomial:  Cusped  Lnd  -  Variation  of 
Independent  Polynomials 
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Figure  10  -  "Quadratic"  Polynomial:  Cusped  End  -  Inflection 

Point  Condition 


Figure  11  -  "Square  Root"  Polynomial:  Rounded  End  -  Variation  of 

Independent  Polynomials 


Figure  12  -  "Square  Root"  Polynomial:  Rounded  End  -  Permissible  Range  of 

Parameters  2r  and  kj 
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Figure  14  -  Ordinary  Polynomial:  Pointed  End  -  Variation  of 
Independent  Polynomials 
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Figure  15  -  Ordinary  Polynomial:  Pointed  F.nd  -  Permissible  Range 

of  Parameters  s  and 


Figure  16  -  Ordinary  Polynomial:  Pointed  Fnd  -  Inflection  Point  Condition 
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Figure  17  -  Ordinary  Polynomial:  Cusped  End  -  Variation  of 
Independent  Polynomials 


Figure  18  -  Ordinary  Polynomial:  Cusped  End  -  Permissible  Rang 


Parameters  k  and  k. 


Figure  19  -  Ordinary  Polynomial:  Cusped  Lnd  -  Inflection  foint  Condition 
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Figure  20  -  "Cubic"  Polynomial:  Flat  Face  -  Variation  of 
Independent  Parameters 
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Figure  21  -  "Cubic"  Polynomial:  Flat  Face  -  Permissible  Range 

of  Parameter  k. 
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